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Abstract 

In this paper, we consider a general cooperative wireless sensor network (WSN) with multiple hops and the problem of 
channel estimation. Two matrix-based set-membership algorithms are developed for the estimation of the complex matrix channel 
parameters. The main goal is to reduce the computational complexity significantly as compared with existing channel estimators 
and extend the lifetime of the WSN by reducing its power consumption. The first proposed algorithm is the set-membership 
normalized least mean squares (SM-NLMS) algorithm. The second is the set-membership recursive least squares (RLS) algorithm 
called BEACON. Then, we present and incorporate an error bound function into the two channel estimation methods which 
can adjust the error bound automatically with the update of the channel estimates. Steady-state analysis in the output mean- 
squared error (MSE) are presented and closed-form formulae for the excess MSE and the probability of update in each recursion 
are provided. Computer simulations show good performance of our proposed algorithms in terms of convergence speed, steady 
state mean square error and bit error rate (BER) and demonstrate reduced complexity and robustness against the time-varying 
environments and different signal-to-noise ratio (SNR) values. 
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I. Introduction 



Recently, there has been a growing research interest in wireless sensor networks (WSNs) because their unique features allow 
a wide range of applications in the areas of military, environment, health and home |1|. They are usually composed of a large 
number of densely deployed sensing devices which can transmit their data to the desired user through multihop relays ftZ\. 
Low complexity and high energy-efficiency are the most important design characteristics of communication protocols [0 and 
| physical layer techniques employed for WSNs. The performance and capacity of WSNs can be significantly enhanced through 
■^j- exploitation of spatial diversity with cooperation between the nodes J2j- In a cooperative WSN, nodes relay signals to each other 
OC in order to propagate redundant copies of the same signals to the destination nodes. Among the existing relaying schemes, the 
C*") ■ amplify-and-forward (AF) and the decode-and-forward (DF) are the most popular approaches @]. Due to limitations in sensor 
(f) | node power, computational capacity and memory (TJ, some power-constrained relay strategies ||5], 10 and power allocation 
■ methods Q have been proposed for WSNs to obtain the best possible SNR or best possible quality of service (QoS) at the 
^""^ destinations. Most of these ideas are based on the assumption of perfect synchronization and available channel state information 
. . ■ (CSI) at each node |lj. Therefore, more accurate estimates of the CSI will bring about better performance in WSNs. 
^ | The normalized least mean squares (NLMS) estimation method is appropriate for WSNs due to its simplicity. However, 
the main problem of the NLMS is that the tradeoff between convergence speed and steady state performance is achieved 
through the introduction of a step size ifTTI . It is not possible to achieve the best solution on these two aspects using a 
C\3 conventional NLMS estimation method. Channel estimation with the NLMS algorithm can be improved by introducing the 
set-membership filtering (SMF) framework lfl2l which modifies the objective function of the NLMS algorithm. It specifies an 
error bound on the magnitude of the estimation error, which can make the step size adaptive. Therefore the SM-NLMS channel 
estimation method can achieve good convergence and tracking performance for each update. A SM-NLMS channel estimation 
algorithm for cooperative WSNs is proposed in IIPJl . Compared with the NLMS channel estimation method, the RLS channel 
estimator can provide better performance in terms of the convergence speed and steady state |[TT1 . However, it is not suitable 
for WSNs due to its high computational complexity IfTTI . In order to overcome this shortcoming, the SMF framework can 
be also introduced to devise a computationally efficient version of the conventional RLS channel estimation method, called 
BEACON channel estimation. It can be considered as a constrained optimization problem where the objective function is the 
least squares (LS) cost function and the constraint is a bound on the magnitude of the estimation error. As a result, an adaptive 
forgetting factor can be derived to achieve the optimal performance for each update. Most importantly, the set-membership 
(SM) algorithms possess a feature that allows updating for only a small fraction of the time, expressed as the update rate (UR). 
Therefore, the UR of the two SM channel estimation algorithms decreases due to the data-selective update which can reduce 
the computational complexity significantly and extend the lifetime of the WSN by reducing its power consumption. 



The biggest issue for the SM channel estimation is the appropriate selection of the error bound, because it has a critical effect 
on the estimation performance. For SM-NLMS channel estimation, the extreme settings of the bound, namely, overbounding 
(the error bound being too large) and underbounding (the error bound being too small) will result in performance degradation 
lfT4l . IfTSl . In practice, the bound depends on the environmental parameters such as the SNR. It is very difficult to determine the 
optimal error bound accurately because there is usually insufficient knowledge about the underlying system. For the BEACON 
channel estimation, the value of the error bound can be varied to trade off achievable performance against computational 
complexity (17). A higher error bound would result in lower UR but worse performance. For WSNs the aim is to achieve 
an acceptable CSI quickly with low power consumption. Therefore, the bound for BEACON channel estimation should be 
adjusted to ensure good estimation performance, lower computational complexity and a low UR. Also, the required error bound 
may be time variant due to changing environmental conditions. 

In this paper, we develop two matrix-based SM algorithms for channel estimation in cooperative WSNs using the AF 
cooperation protocol. The major novelty in these algorithms presented here is that they are matrix-based SM channel estimation 
algorithms as opposed to vector-based SM techniques for filtering applications lfl6l . ifTTl . Ifl8l . Therefore we specify a bound 
on the norm of the estimation error vector instead of the magnitude of the scalar estimation error. Then, a novel error bound 
function is introduced to change the error bound automatically in order to obtain optimal performance with the proposed SM 
channel estimation. Furthermore, we propose analytical expressions of the steady-state output excess mean-square error (MSE) 
of the two SM channel estimation methods. Further novelty in this analysis is that we employ the chi-square distribution to 
describe the probability of the update for estimating the channel matrix as opposed to the Gaussian distribution for estimating 
the filter vector |[T9l . fl20l . A key contribution of this paper is the consideration of techniques to reduce the complexity of the 
channel estimation for WSNs. 

This paper is organized as follows. Section II describes the general cooperative WSN system model and its constrained form. 
Section III introduces two conventional channel estimation methods for reference. Section IV proposes two channel estimation 
methods using the SMF framework and presents an error bound function which tunes the error bound automatically. Section 
V contains the analysis of the steady-state output excess MSE and the computational complexity. Section VI presents and 
discusses the simulation results, while Section VII provides some concluding remarks. 



Consider a general m-hop wireless sensor network (WSN) with multiple parallel relay nodes for each hop, as shown in Fig. 
1. The WSN consists of N s sources, Nd destinations and N r relays which are separated into in — 1 groups: N r m,N r (2\, ... 
>N r r m -i)- All these nodes are assumed to be within communication range. We will concentrate on a time division scheme with 
perfect synchronization, for which all signals are transmitted and received in separate time slots. The sources first broadcast 
the N s x 1 signal vector s to the destinations and all groups of relays. We consider an amplify-and-forward (AF) cooperation 
protocol in this paper. Each group of relays receives the signal from the sources and previous groups of relays, amplifies 
and rebroadcasts them to the next groups of relays and the destinations. In practice, we need to consider the constraints on 
the transmission policy. For example, each transmitting node would transmit during only one phase. In our WSN system, we 
assume that each group of relays transmits the signal to the nearest group of relays and the destinations directly. We can use 
a block diagram to indicate the cooperative WSN system with these transmission constraints as shown in Fig. 2. 

Let H s j,^) denotes the N r u^ x N s channel matrix between the sources and the ith group of relays, H,,^) d denotes the 
Nd x N r (i\ channel matrix between the ith group of relays and destinations, and H r (i-i), r (i) denotes the N r ^ x N r fi_i^ 
channel matrix between two groups of relays. The received signal at the ith group of relays (Xj) and destinations (d) for each 
phase can be expressed as: 



II. Cooperative WSN System Model 



Phase 1: 



x l — H s ,r(l)S + V r (i) 

d 1 = U s .dS + vl 



(1) 



(2) 



Phase 2: 



X 2 — H r ( 1 ) ;r ( 2 )AiXi + V r ( 2 ) 

d 2 =H r(1))d A 1 x 1 +\ 2 d 



(3) 



(4) 



Phase i: (i = 2, 3, m — 1) 



X,; — H r ( i _ 1 ) jr ( i )A i _iX i _i + V r (i) 

d* = H^j^^Aj-iXj-i + 



(5) 



(6) 



Cooperative 
Relays 



Destinations 

Nd 




Fig. 1. An m-hop cooperative WSN with N s sources, iVj destinations and N r relays. 




Fig. 2. Block diagram of the cooperative WSN system with transmission constraints. 



Phase m: 



(m-l),dA m _lX m „i + V™ 



(7) 



where v is a zero-mean circularly symmetric complex additive white Gaussian noise (AWGN) vector with covariance matrix 
a 2 I. Aj is a diagonal matrix whose elements represent the amplification coefficient of each relay of the ith group. The vectors 
d 4 and denote the received signal and noise at the destination nodes during the ith phase, respectively. At the destination 
nodes, the received signal can be expressed as: 

d = H d Ay + v d (8) 



where, 
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Here, we use dashed lines to separate the vectors d, \ d and y in order to distinguish between transmissions to the destinations 
in m different time slots. The matrix consists of all the channels between each group of relays and destinations. The matrix 
A consists of the amplification coefficients of all relays. 

In our transmisstion scheme, all the data packets transmitted from the source nodes and relay nodes contain two parts: a 
preamble part with training sequence symbols and another part with data symbols. Please see Fig. 3. The source nodes transmit 
packets and the relay nodes retransmit those packets that contain the identical training sequence symbols which are known at 
the destination nodes. Therefore, we can make use of them for channel estimation at the destination nodes. After the training 
sequence, the channel estimation algorithm is switched to decision directed mode ETI and the detected data symbols are fed 
to the channel estimator. It can continue to estimate and track the channel. Therefore, the channel variation can be tracked 
after the training phase which can yield better results. Furthermore, this decision directed approach can reduce the length of 
the training sequence which increases the bandwidth efficiency of the WSNs. 



Training Sequence 



Data 



Fig. 3. The structure of the packet transmitted from source nodes and relay nodes 



III. Conventional LS and MMSE Channel Estimation 
Consider a channel estimation problem where the output error is defined as: 

e = r-Hs (12) 

where s (N x 1) is the training sequence symbol vector, H (M x N) is the estimated channel matrix and r (M x 1) is 
the received signal vector at the destination. Conventional channel estimation schemes seek to find the channel matrix H by 
minimizing a cost function which is a suitable objective function of the output error vector e. 



A. The LS Channel Estimator 

The least squares (LS) channel estimation minimizes the weighted sum of the squared norm of the error vector ||e|| 2 which 
can be described as: 



H LS (n) = argmin V A n -'||r(i) - H(n)s(l)f 

H(n) z — ' 

v ' 1=1 



(13) 



where A denotes the forgetting factor. Computing the gradient of the argument and equating it to a zero matrix, we obtain the 
LS channel estimator as given by 11241 : 



5> n -'r(0s*(i) 



1=1 



i=i 



(14) 



where (•) and denote the complex-conjugate (Hermitian) transpose and the inverse respectively. The LS estimator has 
a cubic cost with the number of parameters. A complexity reduction is possible by using a recursive procedure that yields the 
RLS algorithm with quadratic cost. 



B. The MMSE Channel Estimator 

The minimum mean square error (MMSE) channel estimation minimizes the expected value of the squared norm of the 
error vector ||e|| 2 which can be described as: 

Ha/msb = argmin £[||r - Hs|| 2 ] (15) 

H 

After some derivation, the MMSE channel estimator is given by l24l : 

Hmmse = R (S H E[H H H}S + Maty 's^H^H] (16) 

where S and R are the training sequence symbol matrix and received symbol matrix respectively during a training period. The 
MMSE channel estimator requires the full a priori knowledge of the channel correlation matrix and the noise variance cr 2 and 
a cubic cost with the number of parameters. 

IV. Set-Membership Channel Estimation 

In contrast with the two conventional channel estimation methods introduced in section III, set-membership (SM) channel 
estimation specifies an upper bound 7 on the norm of the estimation error vector over a model space of interest which is 
denoted as S, comprising all possible received signal pairs (s,r). The SM criterion corresponds to finding H that satisfies: 

||e(H)|| 2 < 7 2 ,V(s,r)eS (17) 

The set of all possible H that satisfy (17) is referred to as the feasibility set and can be expressed as: 

0= H {HeC MxAr : ||r-Hs|| < 7} (18) 

(s,r)£S 

At time instant n, the constraint set C n is defined as the set of all H(n) that satisfy (17) for the received signal pairs (s(n), r(n)): 

C n = {H(n) GC MxN : \\r(n) - U(n)s(n)\\ < 7} (19) 

The idea behind the SM channel estimation is that if the estimated channel at a time instant lies outside the constraint set C„, 
the estimated channel at the next time instant will lie on the closest boundary of C n . Otherwise, there is no need to compute 
and the power consumption can be significantly reduced. This SM approach makes the estimator adapt only in the direction 
that is necessary. 

A. Proposed SM-NLMS Channel Estimation 

The basic update in the LMS Channel Estimation can be written as: 

H(n + 1) = H(n) + fi(n)e(n)s H (n) (20) 

where e(n) = r(n) — H(n)s(n) denotes the a priori error vector at time instant n, and fi(n) is the time-dependent step size. 
Then we can get a posterior error vector: 

g(n) = r(n) - H(n + l)s(n) (21) 



By substituting (20) into (21), we have: 

g(n) = r(n) - (H(n) + fi(n)e{n)s H (n)) s{n) 

= (r(n) - H(n)s(n)) - /^(»e(n)s ff (n)s(n) (22) 
= e(7i) - ^(n)e(n)s H (n)s(n) 

The constraint set is described as: 

||g(n)|| = ||e(n) - Mn)e(n)s H (n)s(n)|| < 7 (23) 

If ||e(n)|| > 7, then the previous solution lies outside the constraint set. We can choose the constraint value ||g(n)|| equal to 
7 so that the new solution lies on the closest boundary of the constraint set. Therefore: 

||g(n)|| = ||e(n)|| |l - fx(n)s H (n)s(n)| = 7 (24) 

Hence the step size at the nth iteration /j,(n) can be expressed as: 

tin) = ZW7 L \ f 1 ~ IHTTn) ( 25 ) 



s H (n)s(n) \ ||e(n 
Finally, we can write the update equation as: 

H(n + 1) = H(n) + M (n)e(n)s ff (n) (26) 

where, 

H(n) = ( ^(n)s(n) I 1 ~ FRir) ' lf H e ( n )ll > (27) 
^ 0, otherwise. 

Equation (27) shows that the estimated channel matrix updates with a specified step size, only when the norm of the estimation 
error vector is larger than a fixed error bound which we set. Otherwise, the step sizes are zeros which means there is no update 
at these time instants. 

B. Proposed BEACON Channel Estimation 

The proposed BEACON channel estimation method can be considered as the following optimization problem: 

n-l 

minimize £ A(n)"-'||r(0 -H(n)s(Z)|| 2 

1=1 (28) 
subject to ||r(n) - H(n)s(n)|| 2 = -f 2 

To solve this constrained optimization problem, we can modify the LS cost function using the method of Lagrange multipliers 
which yields the following Lagrangian function: 

n-l 

C = Mn) n - 1 HI) - H(n) B (0|| 2 + A(n) [||r(n) - H(n)s(n)|| 2 - 7 2 ] (29) 

where A(n) plays the role of both the Lagrange multiplier and the forgetting factor of the LS cost function. By setting the 
gradient of C with respect to H(n) equal to zero, after some mathematical manipulations (see Appendix), we get the desired 
recursive equation for updating the channel matrix H(n): 

H(n) = H(n - 1) + A(n)e(n)k(n) (30) 

and the recursive equation for updating the intermediate variable matrix P(n): 

P(n) = P(n - 1) - A(n)P(n - l)s(n)k(n) (31) 

where e(n) = r(n) — H(n — l)s(n) denotes the prediction error vector at time instant n. The relationship between k(n) and 
P(n - 1) is 

V( \ = S^(n)P(n-l) 
{ ) l + X{n)s H (n)P(n- l)s(n) k ; 

The error vector is: 

e(n) = r(n) - H(n)s(n) (33) 



TABLE I 

Summary of the BEACON Channel Estimation Algorithm 



Initialize the algorithm by setting 
H(0) = 
P(0) = I 

For each instant of time, n=l, 2, compute 
e(n) = r(n) — H(n — l)s(n) 

m = l 3^(^-0. tf|l«WII>-r. 
[ 0, otherwise. 

where G(n) = s H (n)P(n — l)s(n) 

K W _ l + A(«)G(n) 

H(n) = H(n - 1) + A(n)e(n)k(ra) 
P(n) = P(n - 1) - A(n)P(n - l)s(n)k(n) 



By substituting (30) into (33), we have: 

e(n) = r(n 
= r(n 

= e(n 
= e(n 
= e(n 

= e(n; l + A(n)G(n) 
where G(n) = s H (?i)P(n — l)s(n). The constraint set is described as: 

||e(n)|| = ||e(n) : 



— [H(n - 1) + A(n)e(n)k(n)] s(n) 

— H(n — l)s(n) — A(n)e(n)k(n)s(n) 
inn s H (n)P(n - l)s(n) 

— A(n)e(n) 

1 - 



1 + A(n)s ff (n)P(n - l)s(ra) 
G(n) 



l + A(n)G(n) 
A(n)G(n) 



(34) 



1 + X(n)G(n) 
1 



<7 



(35) 



1 + X(n)G(n) 

If ||e(n)|| > 7, then the previous solution lies outside the constraint set. We can choose the constraint value ||e(n)|| equal to 
7 so that the new solution lies on the closest boundary of the constraint set. Therefore: 

1 



|e(n)|| = ||e(n) 



\l + \(n)G(n)\ 

Hence the optimal forgetting factor at the nth iteration can be expressed as: 



7- 



A(n) = 



1 



\\<n)\\ 



- 1 



G(n) V 7 

Table I shows a summary of the BEACON channel estimation algorithm which will be used for the simulations. 



(36) 



(37) 



C. Time-Varying Bound 

In order to obtain the optimal error bound at each time instant, in this section we introduce an error bound function which can 
adjust the error bound automatically with the update of the channel estimate. A similar bound for the SM filtering techniques 
has been described in lfl2l . For channel estimation, the bound is heuristic and employs the CSI parameter matrix and the noise 
variance that should be related with the estimates of interest. It can be expressed as: 



7(n +!) = (!- /3) 7 (n) + B 'V 'a||H(n)|| s '-ct\ 



(38) 



where f3 is the forgetting factor, a is the tuning parameter, and a 2 is the variance of the noise which is assumed to be known 
at the destinations. This time-varying bound is recursive so that it can be used to avoid too high or low values of ||H(n)|| 2 . 

V. Analysis of the Proposed Algorithms 

A. Steady-State Output MSE Analysis 

In this subsection, we investigate the output MSE in the SM-NLMS and the BEACON channel estimation. The received 
signal at time instant n is given by: 

r(n) = H s(n) + n(n) (39) 



(41) 



where Ho (M x N) is the channel matrix needed to be estimated and n(n) is measurement noise which is assumed here to 
be Gaussian with zero mean and variance a\. Defining the channel estimation error matrix as: 

AH(n) = H H(n) (40) 

we can express the output error vector as: 

e(n) = r(n) - H(n)s(n) 

= r(n) - [H - AH(n)]s(n) 
= r(n) - H s(n) + AH(n)s(n) 
= n(n) + AH(n)s(n) 
Therefore, the output MSE expression can be derived as: 

J(n) = E[\\e(n)\\ 2 } 
= E[e H (n)e(n)] 

= E{[n H (?2)+s H (n)AH H (n)}[n{n) + AH(n)i(n)]} 

= E[\\n(n)\\ 2 } + E[s H (n)AK H (n)AK(n)s(n)} 

= Mai + E{tr[s H (n)AH H (n)AK(n)s{n)}} 

= Ma 2 n + tr{E[s H (n)AH H (n)AH(n)s(n)]} 

where tr(-) denotes the trace of a matrix. The property of the matrix trace tr(XY) = ir(YX) will be used in the following 
derivation. From (42), we can define the output excess MSE as: 

J ex (n) = tr{E[s H (n)AH H (n)AH(n)s(n)]} 
= tr{E[s(n)s H (n)AH H (n)AH(n)]} 



(42) 



(43) 



1) For the SM-NLMS: The update equations for the SM-NLMS channel estimation are given by (26) and (27). In (27) 
s H (n)s(n) is equal to Na 2 , where a 2 s is the variance of the pilot signal. By substituting (27) into (26), we can achieve an 
alternative update equation: 



where 



H(n + 1) = H(n) 
l|eo(n)|| = 



Na 2 



e(n)s H (n) 



||e(n)||, if ||e(n)|| > 7 , 
7, otherwise. 

As a consequence, the update equation of the channel estimation error can be expressed as: 



(44) 
(45) 



AH(n + 1) =AH(n) 



1 



Na 2 



1 



7 



Iko(n) 



e(n)s H (n) 
e(n) 



■s ff (n) 



=AHH--Le(r i ) S fl ( ri ) + -^ 2 
Then, we can use (46) to derive the update equation of the output excess MSE in (43) (see Appendix): 



(46) 



J ex (n + 1) = Ma 2 n + 2 1 E 



1 



Iko(n) 



J ex {n)~2 1 E 



||e(n) 



21 



l|eo(n)| 



■ 7 2 £; 



||e(n) 



|2 1 



l|eo(n)| 



From (45), the three expected values in (47) can be expressed as: 



E 



1 



l|eo(n)|| 



E 



1 



|e(n) 



||e(n)|| >7 



Pup (1 P'Wp ) 

7 



||e(n) 



lko(n)| 



E 



= E[\\t(n)\\\\\t(n)\\> 1 ]P up + -E[Mr 



|e(n) 



||e(n)||< 7 ] (1-P,p) 



l|eo(n)| 



1 



P up + -2£;[||e(n)|| 2 |||e(n)|| < 7 ] (1-P« 



7 



where denotes the conditional expected value and P up stands for the probability of update in each recursion. Let: 

1 



Xx = E 



||e(n) 



||e(n)|| >7 



Yi = E [||e(n)|||||e(n)|| > 7] 



(47) 

(48) 
(49) 
(50) 

(51) 
(52) 



Z!=i5[||e(n)|| 2 |||e(n)|| < 7] (53) 

Equation (47) becomes: 

J ex (n + 1) =M o* + pjX^p + 2(1 - P up )] J ex (n) - 27!^ - 2Zi(l - P up ) + 7 2 ^ P + Zi(l - 
=[27X!P up + 2 - 2P up ] J ex (n) - 2rfYiP up - Z x {\ - P up ) + Ma 2 n + ^P up 

During the steady state, J ex (n + 1) — > J ex (n). Therefore, the steady-state output excess MSE expression of the SM-NLMS 
channel estimation is: 

j ( n \ — 2 l Y i p u P + Z^l - P UP ) - Maj - l 2 Pu P 
ex[U> 2 1 X 1 P up -2P up + l 

2) For the BEACON: According to Table I, we can get the update equation of the channel estimation error for the BEACON 
channel estimation which is very similar to (46): 

AH(„) . AH(„ - 1) - «(■>»' W-'> + .C) sH <">7"-" ,56) 

G(n) ||co(n)|| G{n) 



where, 



. ||6(n)||, if ||e(n)||> 7 , 

1 7, otherwise. v 



Following the same steps described for the SM-NLMS channel estimation in the Appendix, we find that the steady-state output 
excess MSE expression of the BEACON channel estimation has the same style as (55): 

J ex (n) = W, + Z 2 (l-F„ p )-l^- 7 % (5g) 



where, 



Xo = E 



2jX 2 P up - 2P up + 1 
'||e(n)|| >7 



(59) 



\e(n) 

y 2 =£[||e(n)|||||e(n)||>7] (60) 
Z 2 = E[\\e(n)\\ 2 \\\e(n)\\ < 7] (61) 

3) The Probability of Update P up : From (27), we can get the relation about the probability of update of the SM-NLMS 
channel estimation: 

P up = Pr{\\e(n)\\ > 7} = Pr{\\e(n)\\ 2 > 7 2 } (62) 

Similarly, for the BEACON channel estimation we just need to use e(n) instead of e(n). It is easy to see that P up depends 
on the distribution of |e(ri)|| 2 . For the estimated channel matrix Ho with size M x N: 

AI 

||e(n)|| 2 = ^($H 2 [ ei (n)]+3 2 [ e4 (n)]) 

(63) 

al^ SH 2 [ ei (n)] 3 2 [e,(n)] 
2 Q al/2 + all* > 

During the steady state, assuming AH(u) — > 0, the linear relationship between e(n),AH(n) and n(n) in (41) shows that the 
distribution of e(n) is typically Gaussian unless a jamming signal with another distribution is present. Therefore we can get 
that the elements of the error vector e(n) have the same distribution with the elements of the noise vector n(n). Recalling that 
9\[rii(n)] and 3[nj(n)] ~ A/"(0, -f), we can express the distribution of (63) by a chi-square random variable with 2M degree 
of freedom as follows: 



a 



2 



e(n)r~-£Xj M (64) 



Therefore, (62) becomes: 
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TABLE II 

Computational Complexity per Update 



Algorithm 


Multiplication 


Addition 


Division 


NLMS 


2MN + N + min{M, N} 


2AIN + N - 1 


1 


SM-NLMS 


AIN + AI + P up (MN + N + min{M. N}) 


AIN + AI - 1 + P U p{AIN + N) 


2 


RLS 


4N 2 + 2MN + N 


3N 2 + 2AIN - N 


2 


BEACON 


N 2 + AIN + M + N + P up (2N 2 + AIN + N + min{M, N}) 


N 2 + AIN + M - 2 + P up (2JV 2 + AIN - N + 2) 


2 



where F(-) is the chi-square cumulative distribution function (CDF) ||251 defined by: 

Y{l/2) 



W) = TL T] X ' 2) (66) 



By substituting (67) and (68) into (66), we can finally obtain: 

J 



(67) 



In (66) Yl{s,x) is the lower incomplete Gamma function: 

T L (s,x) = I t s_1 e _t di 
Jo 

and Tlx) is the gamma function: 

/>oo 

r(a?) = / f-V** (68) 



J? feO= ! ( L, i _ 1 - (69) 



where 2 denotes the number of degrees of freedom. 



5. Computational Complexity Analysis 

Table II lists the computational complexity per update in terms of the number of multiplications, additions and divisions for 
the SM-NLMS and BEACON algorithms and their competing algorithms. The size of the estimated channel matrix is M x N. 
For our cooperative WSN system model, when is chosen as the estimated channel, we can get: 

M = mN d (70) 

and, 

N = N r + N s (71) 

Because the multiplication dominates the computational complexity of the algorithms, in order to compare the computational 
complexity of our proposed algorithms with their competition algorithms, the number of multiplications versus the size of the 
channel matrix performance for each update is displayed in Fig. 4. For the purpose of illustration, we set M equal to N. It can 
be seen that our proposed SM-NLMS and BEACON channel estimation algorithms have a significant complexity reduction 
compared with the conventional NLMS and RLS channel estimation algorithms. Obviously, a lower P up will cause a lower 
computational complexity. Furthermore, assuming the linear MMSE detectors are used in the destination nodes which require 
cubic complexity, we can get the conclusion that the power used for our proposed channel estimation is only a small fraction 
of the power budget of these nodes. 



VI. Simulations 

In this section, we numerically study the performance of our two proposed SM estimation methods as well as the design 
of the optimal error bound. We consider a 3-hop (m=3) wireless sensor network. The number of sources (N s ), two groups 
of relays (N r m , N r (2)) and destinations (Nj) are 2, 4, 4, 3 respectively. We consider an AF cooperation protocol and the 
amplification coefficient of each relay is set to 1 for the purpose of simplification. We choose as our estimated channel 
because it is the most significant and most complex channel among all channels of the WSN system. The quasi-static fading 
channel (block fading) is considered in our simulations whose elements are Rayleigh random variables and assumed to be 
invariant during the transmission of each packet. Also, in order to test our proposed channel estimation algorithms in a time- 
varying environment, we consider a typical fading channel for wireless communications systems, a Rayleigh fading channel, 
which can be modeled according to Clarke's Model 11261 . According to the transmission scheme introduced in Section II, 
during each phase, the sources and each group of relays transmit the QPSK modulated packets with n p symbols among which 
rit are training symbols and rid are data symbols (Note that n p = n t + rid). n p , n t and rid will be specified in the following 
simulations. The noise at the destinations is modeled as circularly symmetric complex Gaussian random variables with zero 
mean. The SNR is fixed at 10 dB. 




A. MSE performance 

Fig. 5 and Fig. 6 show the channel matrix mean square error (MSE) performance of our proposed SM-NLMS and BEACON 
channel estimation methods for the quasi-static fading channel, and compare them with the conventional NLMS and RLS 
channel estimation algorithms. For the SM-NLMS estimator, we choose five fixed error bounds (7) ranging from 0.3 to 
1.1. It can be seen that increasing the error bound makes the update rate (UR) decrease. It means the update is selective 
which can reduce the computational complexity and power consumption. In the case of an error bound equal to 1.1, the 
UR can fall dramatically to 0.0868. The optimal error bound appears between 0.7 and 0.9. In that situation, the SM-NLMS 
channel estimation achieves the fastest convergence speed and lowest steady states. Otherwise, the performance degrades due 
to overbounding or underbounding. For the BEACON estimator, we choose four fixed error bounds ranging from 0.6 to 0.9. 
Also, the minimum-mean-square error (MMSE) channel estimator which requires the full a priori knowledge of the channel 
correlation matrix and the noise variance is used here for reference. It can be seen that a higher value of 7 results in worse MSE 
performance but a lower UR. In the case of an error bound equal to 0.6, the BEACON algorithm outperforms the conventional 
RLS algorithm (with a forgetting factor of 0.998) in terms of convergence speed and steady state with a slightly reduced UR 
(0.9128). When the error bound is increased to 0.8, although its convergence speed is slower than RLS channel estimation, 
the final MSE is comparable with a much lower UR (0.4356). 

Fig. 7 and Fig. 8 illustrate the performance when we apply the time-varying bound (TVB) into the SM-NLMS and BEACON 
channel estimation. For the SM-NLMS estimator, we transmit packets with 1000 (n p ) symbols among which 100 (rit) are 
training symbols and 900 (rid) are data symbols. We set a to 1.5 and j3 to 0.01. The curve of our proposed algorithm lies on the 
optimal position which is very close to the curve of the SM-NLMS with fixed error bound 0.8. Also, its update rate decreases 
further which is our expectation. For the BEACON estimator, we transmit packets with 2000 (n p ) symbols among which 100 
(nt) are training symbols and 1900 (rid) are data symbols. We set a to 3 and j3 to 0.001. Our proposed algorithm can achieve 
very similar performance to the conventional RLS channel estimation with a substantial reduction in the UR. Therefore, the 
computational complexity is significantly reduced. 

The MSE versus SNR performance of the SM-NLMS and BEACON channel estimation methods are displayed with fixed 
error bounds and the proposed time-varying error bounds in Fig. 9 and Fig. 10. In the cases of fixed error bounds, the MSE 
is lower bounded at different values for different error bounds. For the SM-NLMS estimator, a higher SNR needs a specified 
lower error bound to achieve the optimal MSE performance. When the time-varying error bound is applied, the MSE remains 
very close to the optimal values for all SNRs. For the BEACON estimator, when the SNR is larger than a specified value, 
its MSE will become worse. However, when the time-varying error bound (TVB) is applied, it can be observed that the MSE 
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Fig. 5. MSE performance of the SM-NLMS channel estimation of for quasi-static fading channel compared with the NLMS channel estimation. n p =1000, 
nt=100 and n d =900. 
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Fig. 6. MSE performance of the BEACON channel estimation of for quasi-static fading channel compared with the RLS channel estimation. n p =2000, 
n t =W0 and n d =1900. 
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Fig. 7. MSE performance of the SM-NLMS channel estimation with a time-varying bound for quasi-static fading channel. n p =1000, nt=100 and 71^=900. 
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Fig. 8. MSE performance of the BEACON channel estimation with a time-varying bound for quasi-static fading channel. n p =2000, nt=100 and n ( j=1900. 



keeps on decreasing alone with the increase of the SNR. These two figures show the robustness to the SNR variation of our 
proposed algorithms for the quasi-static fading channel. 
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Fig. 9. SM-NLMS channel estimation MSEs versus SNR for both the fixed bound and time-varying bound for quasi-static fading channel. n p =1000, nt=100 
and n d =900. 
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Fig. 10. BEACON channel estimation MSEs versus SNR for both the fixed bound and time-varying bound for quasi-static fading channel. n p =2000, nt=100 
and n d =1900. 



In order to test our proposed channel estimation algorithms in a time-varying environment, we consider a typical fading 



channel for wireless systems, a Rayleigh fading channel, which can be modeled according to Clarke's Model 11261 . Fig. 11 
and Fig. 12 show the MSE performance of our proposed channel estimation algorithms for the time-varying fading channel 
and three different fading rates (normalized Doppler frequency fdT) are used in the simulations: 10~ 5 , 5 x 10 -5 , and 10~ 4 . 
Because of the requirements of low power consumption and the fact that a fast convergence speed of the proposed algorithms 
might help reducing the need for long training sequences for the WSNs, we focus on the performance of packets with 500 (n p ) 
symbols among which 50 (rit) are training symbols and 450 (rid) are data symbols. For the SM-NLMS estimator, our proposed 
algorithm can achieve better performance than the conventional NLMS algorithms for all the three fading rates. Along with 
the increase of the fading rate, the advantage becomes less pronounced and the update rate becomes higher. For the BEACON 
estimator, our proposed algorithm can achieve very similar performance to the conventional RLS algorithms for all the three 
fading rate. (Note that for the conventional RLS algorithms, when increasing the fading rate, we have to lower the forgetting 
factor to get the optimal performance.) Along with the increasing of the fading rate, the update rate becomes higher. Therefore, 
we can conclude that our proposed channel estimation algorithms can work well for the time-varying fading channel and for 
a range of values of fdT. 
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Fig. 11. MSE performance of the SM-NLMS channel estimation for Rayleigh fading channels compared with the NLMS channel estimation. n p =500, nt=5Q 
and 72^=450. 



B. BER performance 

The MSE performance is very useful to give designers an idea of how well channel estimators perform, whereas bit error rate 
(BER) performance is meaningful in practice. Therefore, in this subsection we focus on the BER performance of our proposed 
algorithms. We consider a simulation where the data packets transmitted at the sources nodes have 1000 (n p ) symbols and 
trained with 100 (nt) symbols. Linear MMSE detectors are used in the destination nodes. We choose as our estimated 
channel and other channels are assumed to be known. Quasi-static fading channel are considered. It can be seen from Fig. 1 3 
that our two proposed SM channel estimation algorithms with time varying bound can achieve a similar BER performance to 
their competing algorithms. Also, the BEACON channel estimator has lower BER than the SM-NLMS channel estimator due 
to the higher computational complexity and the use of the second-order statistics. 

C. Verification of the Analysis 

In this subsection, experiments were conducted to validate our analysis of the SM-NLMS and BEACON algorithms. From 
(70) and (71), the two variables M and N using in the section V can be obtained. M-9 and A^=10. First of all, the analysis of 
the probability of update is verified using (65). It can be seen from Fig. 14 that the P up in simulations of the SM-NLMS and 
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Fig. 12. MSE performance of the BEACON channel estimation for Rayleigh fading channels compared with the RLS channel estimation. n p =500, nt=5Q 
and 71^=450. 
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Fig. 13. BER performance of the proposed channel estimation algorithms. n p =1000, nt=100 and 71^=900. 



BEACON channel estimation is close to and lower bounded by the P up from our analysis. The gap between the analytical curve 



and the simulations of two SM channel estimation is due to the approximation made in the analysis. In section V, we assume 
that the channel matrix error AH approaches zero during the steady-state. However, for the SM algorithms it is not accurate 
because the bound set for the output estimation error would enlarge the AH. During the steady-state, the SM-NLMS channel 
estimation has a larger AH than the BEACON channel estimation which therefore causes a larger gap between the analytical 
curve and the simulation. After that we continue to verify the analysis of the steady-state output excess MSE using (55) and 
(58). Because it is difficult to obtain the full-analytical expressions of the conditional expected values X\ , Y\ , Z\ , X2, Yi , Z2, a 
semi-analytical method is used here. It means that the data from the simulations is used to calculate these conditional expected 
values in (55) and (58). In order to lower the effect of the difference between the analytical P up and the simulation P up of the 
SM-NLMS channel estimation, l.lc^ is chosen approximately to take the place of ct 2 in (65) which would produce a more 
accurate AH and P up for the SM-NLMS channel estimation. Fig. 15 and Fig. 16 show the steady-state output excess MSE 
versus 7 2 /(mA r £ ;cr^) of the two channel estimation algorithms. From the figures, it can be seen that the semi-analytical curves 
can match the simulation curves well. Therefore, it can be stated that our analysis is able to predict accurately the output 
steady-state excess MSE for different choices of bound 7. 
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Fig. 14. Analysis of the probability of the update Pup- 



Nil. Conclusions 

Two SM channel estimation methods have been proposed based on time-varying bound for cooperative wireless sensor 
networks. It has been shown that our proposed methods can achieve better or similar performance to conventional NLMS 
and RLS channel estimation, offering reduced computational complexity. Analyses of the steady-state MSE and computational 
complexity are presented for the two channel estimation and closed-form expressions of the excess MSE and the probability 
of update are provided. Furthermore, the incorporation of the time-varying bound function makes it robust to changes in the 
environment. These features are desirable for WSNs and bring about a significant reduction in energy consumption. 

Appendix 

part of the derivations about the proposed beacon channel estimation algorithm 
By setting the gradient of C in (29) with respect to H(n) equal to zero, we have 

= 2J2 Kn) n - 1 [r(Z) - H(n>(0] [s 11 (/)] + 2A(n)[||r(n) - H(n)s(n)]{~s H (n)} = (72) 
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Fig. 15. Steady-state excess MSE analysis for the SM-NLMS channel estimation. 
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Fig. 16. Steady-state excess MSE analysis for the BEACON channel estimation.. 



Therefore, 



H(n) 



n-l 



^2 \(n) n - l s(l)s H (I) + \(n)s(n)s H (n) 



1=1 



= Hn) n ~ l r(l> H (l) + A(n)r(n)s ff (n) 



(73) 



z=i 



Then we can get 



H(n) 



n-l 



53 \(n) n - l r(l)s H (I) + \(i2)r(n)s H {n) 



1=1 



n-l 



53 X(n) n - l s(l)s H (l) + X(n)s(n)s H (n) 



i=i 



Let: 



and, 



4>{n) = 53 A(n)"-'s(Z)s ff (0 + \{n)s(n)s H (n) 



i=i 



Z(n) = 53 A(n) n -'r(0s ff (0 + A(ri)r(n)s ff (n) 



z=i 



Equation (74) becomes: 



H(n) = Z(n)0 _1 (n) 



(74) 



(75) 



(76) 



(77) 



Isolating the term corresponding to I = n — 1 from the rest of the summation on the right-hand side of (75), we may write: 

"n-2 

(78) 



0(n) = 



53 A(7i) n -'s(0s H (Z) + A(n)s( ? i - l)s ff (7i - 1) 



z=i 



+ A(n)s(n)s ff (n) 



The expression inside the brackets on the right-hand side of (78) equals <p{n — 1) assuming the forgetting factor of the cost 
function is close to 1. Hence, we have the following recursion for updating the value of 4>(n): 



4>(n) = 4>{n - 1) + X(n)s(n)s H (n) 

Similarly, we may use (76) to derive the following recursion for updating Z(n): 

Z(n) = Z(n - 1) + A(n)r(n)s H (n) 
Then, using the matrix inversion lemma ifTTI . we obtain the following recursive equation for the inverse of 4>(n): 

A-ifn) = A -U n _,)_ Kn)^\n - l) S (n)s H (n^n)^ 1 (n - 1) 
V 1 ' V K ' 1 + A(n)s»(n)</.- 1 (n- l)s(n) 



For convenience of computation, let: 
and, 

Therefore, we may rewrite (77) and (81) as: 



k(n) = 



P(n) = _1 (n) 

S g (n)P(n- 1) 
l + A(7i)s H (n)P(n- l)s(ra) 



H(n) = Z(n)P(n) 
P(n) = P(n - 1) - A(n)P(ra - l)s(n)k(n) 
Then we substitute (80) and (85) into (84) to obtain a recursive equation for updating the channel matrix H(n): 

H(n) = H(n — 1) — A(n)H(n — l)s(n)k(n ) + A(n)r(n)s ff (n)P(rc) 
By rearranging (83) , we can get: 

k(n) = s H (n)P(n - 1) - A(n)s H (7i)P(n - l)s(n)k(n) 
= s H {n) [P(n - 1) - A(7i)P(n - l)s(n)k(n)] 
= s H (n)P(n) 

Using (87) above, we get the desired recursive equation for updating the channel matrix H(n): 

H(n) = H(n - 1) - A(n)H(n - l)s(n)k(n) + A(n)r(n)k(n) 
= H(n - 1) + A(n) [r(n) - H(n - l)s(n)] k(n) 
= H(n - 1) + A(n)e(n)k(n) 

where e(n) = r(n) — H(n — l)s(n) denotes the prediction error vector at time instant n. II 
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PART OF THE ANALYSIS OF THE PROPOSED SM-NLMS CHANNEL ESTIMATION ALGORITHM 

From (46), the update equation of the channel estimation error is: 



AH(n + 1) = AH(n) 
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No* ||e (n) 
Equation (89) becomes: 

AH(ra + l) = A + B 

From (43), we can get the output excess MSE at time instant n + 1: 

J ea; (7i + 1) = <r{£[s(n + \)% H {n + l)AH H (n + l)AH(n + 1)]} 
= tr{E[s(n)s H (n)AH H (n + l)AH(n + 1)]} 

= -01 + -02 + V'3 + ^4 

Then we analyze each term separately: 
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Finally, we can obtain the update equation of the output excess MSE: 
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